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SIMULATION  OF  A  URANIUM  VAPOR  RELEASE 
IN  THE  HIGH  ALTITUDE  ATMOSPHERE 


I .  Introduction: 

High  altitude  nuclear  explosions  can  lead  to  greatly  enhanced  infrared 
radiation  in  the  earth's  atmosphere  and  can  thereby  have  a  deleterious  effect 
on  defensive  infrared  systems.  This  enhanced  infrared  radiation  has  a  number 
of  sources:  direct  plasma  radiation  from  the  nuclear  burst  created  plasma; 

creation  of  increased  densities  and/or  increased  temperatures  of  natural 
atmospheric  infrared  radiators  such  as  hydroxyl,  nitrogen  oxides  and  carbon 
oxides;  and  also  infrared  radiation  of  weapon  and  vehicle  debris  species  which 
are  vaporized  in  the  explosion.  Recently  it  has  been  recognized  that  urnium 
and  its  oxides  are  likely  to  be  the  dominant  infrared  active  species  at  1  i 
minutes  and  later  following  a  high  altitude  nuclear  detonation  ( Arms t rone , 
1981). 

The  Defense  Nuclear  Agency  has  an  interest  in  developing  methods  to 
predict  the  infrared  radiation  after  a  high  altitude  nuclear  event.  Armstrong 
points  out  those  areas  where  there  Is  a  severe  lack  of  knowledge  required  to 
make  accurate  predictions.  In  order  to  measure  the  spectral  lines  and  the 
strength  of  long  wavelength  infrared  emission  it  is  expected  that  fairly  large 
releases  of  uranium  vapor  will  be  necessary  in  the  earth's  high  atmosphere. 
Only  after  such  measurements  are  made  will  accurate  predictions  become 
possible.  With  this  requirement  in  mind,  we  have  considered  the  scenario 
originally  described  by  Reidy  (1980)  involving  the  prompt  deposition  of  2  kg 
of  uranium  vapor  at  an  altitude  of  200  km.  Current  experimental  release 
techniques  are  short  of  being  capable  of  this  rapid  a  deposition  by  orders  of 
magnitude. 

In  this  memorandum,  we  present  results  for  the  temporal  evolution  of  a 
uranium  vapor  release  at  approximately  200  kilometer  altitude  in  the  earth's 
atmosphere.  The  results  are  obtained  by  numerical  simulations  using  a  simple 
model  and  can  be  used  to  plan  such  an  experimental  release.  The  results 
describe  both  the  chemical  and  hydrodynamic  evolution  for  one  experimental 
scenario  for  time  periods  of  a  few  tens  of  seconds.  We  predict  a  uranium  ion 
cloud  extended  along  the  field  but  narrow  perpendicular  to  the  earth's 
magnetic  field.  The  Ion  densities  are  high  for  minutes  after  release  with  the 
Manuscript  approved  August  18,  1988. 
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II .  Simulation  Model: 

We  consider  a  simple  one-dimensional  model  of  the  uranium  release.  The 
temporal  evolution  of  each  species,  denoted  a,  is  described  by  equation  1. 

3n  3n  3n  3n 

— —  -  — -  hydro  +  — —  d if  f  +  - — —  chem  (1) 

3t  3t  3t  3t 

The  first  term  on  the  right  is  the  change  in  the  concentration  from 
hydrodynamic  transport.  The  second  term  describes  diffusion  of  ions  through 
ions  and  neutrals  through  neutrals.  The  third  term  treats  chemical 
transformation  of  the  individual  species. 

The  hydrodynamic  transport  of  ions  and  neutrals  is  done  separately.  The 
transport  equations  are: 


continuity 
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Here  s  denotes  either  the  ions  or  the  neutrals  and  t  denotes  the  other.  The 
symbols  p  ,v_,  p,  T,  and  q  are  density,  flow  velocity,  pressure,  temperature, 


and  artificial  viscosity  respectively.  y  is  the  ratio  of  specific  heats  and 
a  the  collisional  coefficient  of  friction. 

In  this  simple  one  dimensional  model  of  transport  we  treat  two  specific 
cases;  transport  perpendicular  and  transport  parallel  to  the  ambient 
geomagnetic  field.  Perpendicular  to  B  the  neutrals  are  transported 
spherically  and  the  ions  are  frozen  to  the  unmoving  field.  Parallel  to  B,  the 
neutrals  are  transported  spherically  and  the  ions  are  transported  along  a 
cartesian  axis  parallel  to  the  field.  The  two  cases  of  the  simple  one¬ 
dimensional  model  should  give  reasonable  estimates  for  the  size  of  a  release 
cloud  in  the  atmosphere.  Since  the  transport  model  allows  for  diffusion 
between  ions  and  neutrals  our  diffusion  model  only  needs  to  consider  neutrals 
and  ions  diffusing  through  themselves. 

The  diffusion  equation  in  a  multi  constituent  gas  can  be  expressed 
(Burgers  1969) 


a  "  V  =  T(VP«"~  Vp)' 

A 


(5) 


where  p  is  the  total  density.  An  exact  solution  for  diffusion  can  be  obtained 

by  inversion  of  (5)  and  substitution  of  the  species  velocities  into  a  set  of 

coupled  continuity  equations  for  each  specie.  For  our  model  here  however  we 

approximate  (5)  and  obtain  a  decoupled  set  of  continuity  equations.  We  assume 

a  =  a  ,  a  constant,  which  allows  us  to  express  the  effects  of  diffusion  as 
010 


at 


=  -v» ( K V  p  /p), 

01 


(6) 


where  <  is  a  diffusion  coefficient  depending  on  o  .  We  have  tested  this 
simplified  diffusion  model  and  have  obtained  reasonable  results.  The 
simplification  of  diffusion,  while  not  totally  accurate,  yields  results  which 
are  within  a  factor  of  two. 


The  effects  of  chemistry  on  the  species  are  computed  in  a  series  of  rate 
equations 
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The  reactions  of  concern  and  their  rate  coefficients,  kgt, 
1  which  are  taken  from  Armstrong  (1981)  and  Archer  (1982). 


are  shown  in  Table 


The  final  part  of  the  model  is  the  specification  of  the  neutral 
atmosphere.  We  have  taken  species  concentrations  and  a  temperature 
appropriate  to  moderate  solar  activity,  daytime,  200  km  altitude  conditions; 
they  are  [0]  =  4. 4  x  10^  cm-'*,  [^1  =  A. 2  x  10^  cm-*,  and  [07]  = 

2.6  x  10®cm-*.  The  neutral  atmospheric  temperature  is  0.1  ev. 

III.  Results: 

We  have  solved  the  model  equations  1,2, 3, 4, 6,  and  7  for  a  uranium  vapor 

release  at  200  km  altitude  in  the  atmosphere.  We  have  assumed  a  prompt 

2  2 

release  of  2  kg  of  atomic  uranium,  approximately  5.  x  10  atoms.  The 
temperature  of  the  uranium  vapor  is  taken  as  0.276  ev.  Three  sets  of 
solutions  were  obtained  with  the  diffusion  coefficient,  k,  as  a  parameter  to 
study  the  effects  of  diffusion  on  the  release  evolution.  In  the  first,  or 

2 

nominal,  case  <  was  taken  as  0.1  km  /sec,  a  nominal  diffusion  coefficient  at 
200  km  altitude.  For  an  appropriate  gradient  scale  length  of  100  m  the 
diffusion  velocity  is  100  m/sec.  This  value  is  about  17%  of  the  sound  speed 
at  that  altitude.  We  have  also  obtained  results  for  diffusion  coefficients 
1/3  and  3  times  the  nominal  value.  The  slower  diffusion  rates  lead  to  very 
slow  oxidation  of  the  uranium  vapor  and  are  therefore  uninteresting.  In  the 
following  sections  we  show  results  for  nominal  diffusion  and  faster  diffusion. 

Figures  1  through  14  show  results  for  the  simulation  perpendicular  to  the 
geomagnetic  field  with  nominal  diffusion.  Figure  1  shows  the  initial 
conditions  of  the  release.  The  vertical  axis  plots  the  logarithm  to  the  base 
10  of  the  gas  concentration  while  the  horizontal  axis  measures  radius  of  the 
release  in  km.  The  uranium  vapor  has  a  radius  of  200  m  and  has  displaced  the 
ambient  atmospheric  species.  We  have  tested  a  number  of  initializations  and 
find  that  the  results  after  a  few  seconds  are  insensitive  to  details  of  the 
Initial  conditions. 

Figures  2,  3,  and  4  show  the  neutral  gas  concentration  at  2,  5,  and  10 
seconds  respectively.  Initially  both  atomic  and  molecular  oxygen  are  burned 
out  of  the  release  cloud  in  oxidizing  the  uranium  gas.  The  atomic  oxygen  is 
able  to  diffuse  into  the  cloud  by  about  5  seconds  but  the  molecular  oxygen 
continues  to  be  consumed  as  it  diffuses  in  for  times  longer  than  10  seconds. 


Table  l 


Chemical  Reaction  Rare  Coefficients 


(1) 

U  +  0  •»  U0f  t-  e~ 

5  x 

10  ^  cm”* / sec 

(2) 

U  a-  0  *  U0+  +  e” 

3.8 

x  10  1  ~  cinVsec 

(3) 

n  +  o9  »  UO  +  0 

1.8 

x  10"1'’  cmVsec 

(4) 

UO  +  0  +  UOt  1-  e 

4  x 

in-10  cm^soe 

(5) 

UO  +  09  +  U0.j+  0 

2  x 

10-0  cm1.' sec 

(6) 

uof+  02>  UO*  +  0 

4  x 

-n  8 , 

10  cm  -sec 
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The  different  behaviour  of  atomic  and  molecular  oxygen  is  due  primarily  to 
their  different  concentrations  in  the  ambient  atmosphere;  there  just  is  not 
enough  molecular  oxygen  around  to  rapidly  fully  oxidize  the  cloud. 

The  resulting  ion  cloud  concentrations  are  presented  in  Figures  5  through 
1 1  for  post  release  times  of  2  to  50  seconds.  The  U0+  component  forms  in  a 
shell  like  structure  at  a  radius  of  about  200  meters  from  the  release  point. 
The  shell  is  formed  at  a  burning  front  between  outward  diffusing  uranium  vapor 
and  inward  diffusing  atomic  oxygen,  and  the  ions  are  frozen  in  the  geomagnetic 
field  at  their  point  of  oxidation.  The  initial  oxidation  of  uranium  gas  to 
U0+  is  complete  by  5  seconds.  Further  burning  of  U0+  to  U0*  by  molecular 
oxygen  then  takes  place  in  an  outer  shell  which  slowly  moves  inward  as 
unconsumed  molecular  oxygen  is  able  to  diffuse  inward.  The  complete  oxidation 
of  uranium  vapor  to  UO*  takes  between  40  and  50  seconds.  The  final  ion 
densities  of  10  cm  in  the  shell  are  higher  than  would  occur  in  an  actual 
release  because  this  one-dimensional  model  does  not  allow  for  diffusion 
parallel  to  the  magnetic  field.  From  studying  our  one  dimensional  parallel 
release  results  we  would  expect  the  maximum  ion  cloud  densities  to  be  about 

q 

10  cm  for  this  release. 

Figures  12,  13,  and  14  show  the  results  for  U0+,  U0*,  and  0., 
concentrations  respectively  in  a  different  graphical  format.  We  have  plotted 
contours  of  concentration  versus  radius  and  time  after  release.  Here  the 
shell  structure  of  the  ion  cloud  is  clearly  seen.  The  U0+  forms  a  shell  at 
200  meters  radius  immediately  after  release.  The  U0*  shell  initially  forms 
at  about  400  meters  and  gradually  builds  inward  as  oxidation  takes  place.  We 
are  also  able  to  see  the  cavity  formed  in  the  ambient  which  is  only  slowly 
filled  by  diffusion  from  the  outside. 

K  similar  set  of  results  for  the  evolution  of  the  uranium  release 
parallel  to  the  geomagnetic  field  is  shown  in  Figures  15  through  24.  The 
initial  conditions,  which  are  identical  to  the  previous  case  are  shown  in 
Figure  15;  note  the  change  in  scale  on  the  horizontal  axis.  The  neutral 
concentrations  for  2,  5,  and  10  seconds  after  release  are  shown  in  Figures  16, 
17,  and  18  respectively.  The  ion  concentrations  for  the  same  time  intervals 
after  release  are  shown  in  Figures  19,  20,  and  21.  The  qualitative  nature  of 
the  temporal  and  spatial  evolution  is  like  that  discussed  previously  for  the 
perpendicular  release.  The  primary  difference  in  the  results  is  that  here  the 


ion  cloud  is  able  to  readily  diffuse  outward.  This  parallel  diffusion  allows 
the  U0+  to  become  very  extended  along  the  geomagnetic  field.  The  L'0+  is 
however  still  confined  to  the  release  region  as  it  is  consumed  by  oxidation  by 
molecular  oxygen  as  it  attempts  to  diffuse  parallel  to  the  magnetic  field. 

Figures  22,23  and  24  show  contours  of  U0+,  UO^  ,  and  0^  concentration 
respectively  for  the  parallel  simulation.  Again  the  results  are  similar  Lo 
those  for  the  perpendicular  simulation.  The  formation  of  an  ion  cloud  shell 
about  the  release  point  is  again  clearly  seen.  The  shell  is  broader  in  this 
case,  and  for  U0+  diffusion  outward  is  very  noticeable.  The  molecular 
oxygen  is  depleted  in  the  release  location,  and  fot  tills  case  the  burned  out 
region  is  very  much  larger  as  the  Lons  diffuse  outwird. 

The  results  for  the  uranium  release  simulations  with  3  times  faster 
diffusion  are  shown  in  Figures  25  through  2b.  The  initial  conditions  are  the 
same  as  used  for  normal  diffusion.  Figures  23  and  2h  present  results  of  the 
perpendicular  simulation  10  seconds  after  release  for  the  neutrals  and  ions 
respectively.  Figures  27  and  28  show  results  for  the  parallel  simulation. 
The  results  in  both  instances  are  qualitatively  similar  to  those  represented 
previously.  The  major  quantitative  differences  are  a  more  apid  oxidation  of 
U0+  to  U0+  allowed  by  faster  diffusion  of  0o  into  the  ion  cloud,  and  for  the 
paralle1  case,  a  more  rapid  diffusion  of  U0+  along  the  geomagnetic  field. 
Nevertheless  the  time  scale  for  the  complete  oxidation  of  U  to  U<)_  is  still 
greater  than  10  secs. 

IV.  Cone  Ins  ions 

The  results  of  the  numerical  simulation  of  a  uranium  vapor  release  at 
about  200  km  altitude  have  demonstrated  a  very  consistent  and  understandable 
morphology.  For  the  2  kg  prompt  release  considered  here,  it  is  clearly 
demonstrated  that  the  uranium  vapor  will  completely  burn  out  the  molecular 
oxygen  of  the  ambient  neutral  atmosphere.  Moreover  diffusion  of  molecular 
oxygen  from  outside  into  the  release  cloud  is  not  sufficient  to  rapidly 
oxidize  the  cloud  to  its  final  state,  UO^  .  We  have  shown  that  the  time 
scale  for  complete  oxidation  Is  a  few  tens  of  seconds. 

It  might  be  thought  that  increasing  the  altitude  of  release  would 
alleviate  the  slow  oxidation  of  the  cloud  owing  to  faster  diffusion  rates. 
Unf ortunately  this  is  probably  not  the  case  since  the  atmospheric  molecular 
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oxygen  density  decreases  rapidly  with  altitude  and  the  burn  out  would  be  more 
severe.  It  would  appear  that  the  only  solution  t<  the  slow  oxidation  problem 
would  be  to  release  the  uranium  vapor  gradually.  However,  in  this  case  one 
would  not  have  a  compact  cloud  for  measurement  but  would  have  a  long  thin 
trail  of  oxidized  uranium  ions  across  the  sky. 

The  morphology  of  the  release  simulated  here  has  a  clearly  defined  shell 
structure.  The  ion  cloud  is  also  greatly  extended  along  the  geomagnetic 
field.  At  early  times  after  release  the  Ions  form  a  double  shell  with  UO* 
ions  forming  outside  the  U0+  shell.  At  later  times  there  is  a  single  shell  of 
UO^  ions  with  a  diameter  of  about  1  km  perpendicular  to  the  magnetic  field 
and  a  length  along  the  field  of  more  than  10  km.  The  simulation  indicates 

Q  Q 

that  the  ion  density  in  the  shell  can  be  expected  to  be  10  to  10  cm  .  Til 
is  a  very  high  density  when  compared  to  the  well  studied  barium  vapor  release 
where  the  Ion  densities  are  normally  2  orders  of  magnitude  lower. 

These  very  high  Ion  densities  can  have  important  implication  fir  gradiei 
drift  structuring  for  the  uranium  release  ion  cloud.  For  barium  releases,  i  K 
has  been  demonstrated  that  ion  clouds  which  have  a  large  conductivity  ratio 
compared  to  the  background  ionosphere  form  striations  very  slowly  (Hinson, 

1975  ).  Therefore,  we  would  expect  that  a  uranium  vapor  release  would  create 
a  slowly  structuring  ion  cloud.  This  situation  may  be  somewhat  alleviated  by 
conducting  the  uranium  release  during  daytime  with  higher  ambient  ionospheric 
conductivities.  However,  on  balance  the  uranium  ion  cloud  should  only  slowly 
form  striations. 

In  summary,  we  conclude  that  a  prompt  uranium  vapor  release  would  result 
In  a  UO^  Ion  cloud  In  a  few  tens  of  seconds.  The  ion  cloud  would  have  a 
shell  like  structure  and  be  compact  perpendicular  to  the  magnetic  field  and 
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very  elongated  along  the  field.  The  Ion  concentrations  would  range  from  10 

9  —3 

to  10  cm  and  would  develop  striated  structure  very  slowly. 
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Figure  10.  Ion  gas  concentrations  at  40  set-onus  after  release. 
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Figure  15.  Initial  conditions  for  simulation  paral 
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Figure  17.  Neutral  gas  concentrations  at  r>  sivumls  aft 
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Figure  18.  Neutral  gas  concentrations  at  10  seconds  after  release 
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Figure  20.  Ion  concentrations  at  ^  seconds  after  release 
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Figure  2f>.  Ion  gas  concentrations  at  10  seconds  For  simulation  perpendicular 
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1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  S7106 
01CY  R.  STKLLINCWERF 
OICY  M.  ALME 
01CY  L.  WRIGHT 

MITRE  CORPORATION,  THE 

P.O.  BOX  208 

BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANSTERN 
OICY  ATTN  C.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATF.  RESEARCH  PARK 

1820  DOLLY  MADISON  BLVD 

MCLEAN,  VA  22101 
OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

P AC  I F I C-S  l  F.RRA  RESEARCH  CORP 

12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY ,  CA  90291 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KAR7.AS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 
OICY  ATTN  L.  DeRAND 
OICY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 
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SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
01CY  ATTN  R.  TSUNODA 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CKESNUT 
01CY  ATTN  CHARLES  L.  RINO 
0  ICY  ATTN  WALTER  JAYE 
01CY  ATTN  J.  VICKREY 
01CY  ATTN  RAY  L.  LEADABRANO 
01CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
01CY  ATTN  R.  LIVINGSTON 
0 ICY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 

01CY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  6  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R.  K.  PLEBUCH 
0 ICY  ATTN  S.  ALTSCHULER 
0  ICY  ATTN  D.  DEE 
01CY  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
0 ICY  ATTN  W.  REIDY 
0 ICY  ATTN  J.  CARPENTER 
0 ICY  ATTN  C.  HUMPHREY 


